Environment Assisted Precision Measurement 



G. Goldstein^, P. Cappellaroi^^.s^ Maze\ J. S. Hodges^'^, L. Jiang^^^^ A. S. S0rensen4, M. D. Lukini-^ 
^Department of Physics, Harvard University, ^ Harvard- Smithsonian Center for Astrophysics, Cambridge MA 02138 USA 
'^Nuclear Science and Engineering Dept., Massachusetts Institute of Technology, Cambridge MA 02139 USA 
^ QUANTOP, Niels Bohr Institute, Copenhagen University, DK 2100, Denmark and 
^Institute of Quantum Information, California Institute of Technology, Pasadena CA 91125 USA 



o 
o 

(N 
o 

Q 



Oh, 

a ■ 



> 

00 

o 
p 

o 
o 



X 



We describe a method to enhance the sensitivity of precision measurements that takes advantage 
of a quantum sensor's environment to amplify its response to weak external perturbations. An 
individual qubit is used to sense the dynamics of surrounding ancillary qubits, which are in turn 
affected by the external field to be measured. The resulting sensitivity enhancement is determined 
by the number of ancillas that are coupled strongly to the sensor qubit; it does not depend on the 
exact values of the coupling strengths and is resilient to many forms of decoherence. The method 
achieves nearly Heisenberg-limited precision measurement, using a novel class of entangled states. 
We discuss specific applications to improve clock sensitivity using trapped ions and magnetic sensing 
based on electronic spins in diamond. 



Precision measurement is among the most important 
applications of resonance methods in physics. For exam- 
ple, quantum control of atomic systems forms the physi- 
cal basis of the world's best clocks. Ideas from quantum 
information science have been used to demonstrate that 
quantum entanglement can enhance these measurements 
At the same time a wide range of quantum systems 
have been recently developed aimed at novel realizations 
of solid state qubits. Potentially such systems can be 
used as quantum measurement devices such as magnetic 
sensors with a unique combination of sensitivity and spa- 
tial resolution [3,0, 12, 3^. In this Letter, we describe a 
novel technique that makes use of the sensor's local en- 
vironment as a resource to amplify its response to weak 
perturbations. We shall use solid state sensors and ion 
clocks as examples. 

The purpose of quantum metrology is to detect a small 
external field, coupled to the sensor by a Hamiltonian: 
= b{t)KSz, where is the spin operator of the quan- 
tum sensor. Here, b(t) can be an external magnetic field 
or the detuning of a laser from a clock transition while k 
is the spin's coupling to the field. The working principle 
of almost any quantum metrology scheme can be reduced 

where the field is mea- 



to a Ramsey experiment 1^, 1 
sured via the induced phase difference between two states 
of the quantum sensor. The figure of merit for quantum 
sensitivity is the smallest field Sbmin that can be read out 
during a total time T. For a single spin 1/2, if the sensing 
time is limited to r (e.g. by environmental decoherence) 
then: dbrmn ^ 

In many cases the external field also acts on the sen- 
sor's environment, which normally only induces decoher- 
ence and limits the sensitivity. Here we show that in 
some cases the environment can instead be used to en- 
hance the sensitivity. For generality we will illustrate the 
key ideas using the central spin model (Fig. la). In this 
model a central spin (which can be prepared in a well de- 
fined initial state, coherently manipulated and read out) 
is coupled to a bath of dark spins that can be polar- 



ized and collectively controlled, but cannot be directly 
detected. The system is described by the Hamiltonian 
H ^ Hb + Hint, with 



Hb = b{t){KS, 



mnt = \l){l\Y.\rIl, (1) 



where are the couplings between sensor and environ- 
ment spins, while n and ^ are couplings to the external 
field of the central and dark spins respectively. Here |0), 
|1) and Sz refer to the central spin while |t), II to the 
dark spins (and we set h = 1). We consider two cases. 
In the first one Hi^t can be turned on and off at will 
and is much larger than any other interaction in the sys- 
tem (e.g. a laser-mediated ion interaction). In the second 
case. Hint is intrinsic to the system and of the same order 
of magnitude as the relevant sensing time (e.g. dipole- 
dipole interactions between solid state spins) . In all cases 
we will assume collective control over the dark spins. 

To illustrate the sensing method we consider first the 
idealized case where the couplings between the central 
and the dark spins can be turned on and off at will and 
the dark spins are initialized in a pure state |tt ••• t)- 
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FIG. 1: Ideal Model, (a) A central spin coupled to a spin 
bath, (b) Ideal measurement procedure. Gates labeled \Ix 
represent controlled rotations e~'-^^= of the dark spins, ob- 
tained via the Hamiltonian Hint- The gates bl^ represent 
rotations e'^^'" due to the external field b. The central spin 
undergoes a spin-echo before measurement. 

Consider the circuit in Fig. 1(b). First, the central spin 
is prepared in an equal superposition of the two internal 
states ^ |0) -f |1). Then Hint is rotated to the x-axis and 
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applied for a time r. This induces controlled rotations of 
the dark spins, resulting in an entangled state 

(|0)|t...t> + |l)l^i-.¥'^»/V2, (2) 

where \Lpi) = cos ((^i) |t) — i sin((pj) ||) with Lpi = XiT. 
This state is then used to sense the magnetic field. We 
assume for now k = and define = ^ /□ dtb{t). Un- 
der the action of the magnetic field the state evolves to 
(|0)|t...) + |l>l^i---^Jv))/\/2 with 

IV'i) = cos(p, It), - isin(^,e2»«<i ||)^ 

where {Lpi\ifj^'^ = 0. The central spin is then flipped by 
a TT-pulse and another control operation with Hint along 
X is applied. If the field (and thus 9d) were zero the 
interaction between central and dark spins would then 
be refocused, corresponding to decoupling the sensor spin 
from the environment as in a spin-echo. For metrology 
purposes, to first order, the effect of a small field is to 
introduce a phase difference $ « 29 d siv? [ipi) between 
the states of the dark spins, depending on the state of 
the central spin (while the terms oc \'-P't) only contribute 
to second order in 9d)- After the sensor spin is rotated 
around the x-axis, this yields an additional contribution 
to the probability of finding the spin in the |l)-state: 
Pi w (H- $ -I- 9s) /2 -I- 0(&^) (where we reintroduced the 
phase acquired by the sensor spin alone, 9^ = n JJ" dt b{t), 
which can be simply added as commutes with the rest 
of the Hamiltonian) . 

While the signal is enhanced by a factor cx the quan- 
tum projection noise remains the same as we still read 
out one spin only. The minimum field that can be mea- 
sured in a total time T is then: 

Sbrn.n = « (3) 

where n is the total number of dark spins. The linear 
scaling in n of the phase $ can be achieved in principle 
for any distribution of A^'s, since we can always choose 
a duration r such that (sin^ (Air)) > ^, leading to or- 
der one contribution from each spin. Thus we are able to 
perform Heisenberg-limited spectroscopy despite the fact 
that the precise form of the entangled state ^ is uncon- 
trolled, and may not even be known to us [20|. This 
considerably relaxes the requirements for entanglement 
enhanced spectroscopy as compared to known strategies 
involving squeezed or GHZ-like states. We next discuss 
two experimental implementations that approximate this 
idealized scheme: quantum clocks with trapped ions and 
spin-based magnetometry. 

To reach high precision in quantum clocks, the ions 
must posses several characteristics: a stable clock transi- 
tion, a cooling cycling transition, good initial state prepa- 
ration and reliable state detection. It is then convenient 



to use two species of ions in the same Paul trap [24|, |30| : 
The spectroscopy ions (e.g. '^'^Al'^) provide the clock 
transition while the logic ion (e.g. ^Be~^) fulfills the 
other requirements. Although inspired by a similar idea 
as the one proposed here, experiments using two ion 
species have been so far limited to just one spectroscopy 
ion ll^, [iil - fsol l ; with our method the number of spec- 
troscopy (dark) ions can be increased. 

Specifically, by using multichromatic gates |2ll| o ne can 
implement the Hamiltonian Hint in Eq. ^ [2o!| (such 
multichromatic gates are known to be much more robust 
to heating noise than Cirac-Zoller gates [2l| that have 
been used so far). We can then use the method pre- 
sented here to transfer the phase difference due to the 
detuning of several spectroscopy ions onto a single logic 
ion, which can be read out by fluorescence. In princi- 
ple, we achieve Heisenberg-limited sensing of the clock 
transition without individual addressability of the spec- 
troscopy ions or producing GHZ states. Importantly, we 
achieve this even if the spectroscopy ions have different 
couplings to the logic ion, which will be the case in a trap 
with different ion species due to the absence of a common 
center of mass mode. 

We next briefly discuss the effects of decoherence on 
this method. It has been argued (see e.g. [3l|), that the 
coherence time reduction for a n-spin entangled state re- 
duces the sensitivity to roughly that of spectroscopy per- 
formed on n individual spin, scaling as ^/rl, as opposed to 
the ideal scaling c>c n derived here (See Ref. for a dif- 
ferent scaling in the case of atomic clocks) . Our scheme 
would obtain an improvement cx -yAi with respect to cur- 
rent experimental realizations where only a single ancil- 
lary ion is available, even in the decoherence model of 
Ref. [sij, where each spin undergoes individual Marko- 
vian dephasing. Furthermore, this decoherence model 
is not so relevant in present setups, as technical noise 
during the gates and imperfect rotations are dominant 
for traps with many ions. Our method is highly robust 
to static repetitive imperfections during the gates, lead- 
ing to further improvement depending on the exact noise 
specifications [26(|. 

In many physical situations short bursts of controlled 
rotations, as used above, are not available. Instead, the 
couplings between the central and dark spins are always 
on and their exact strength is unknown. Examples of 
such systems are solid-state spin systems used for mag- 
netometry 0, 3, m, ll^l- Still, it is possible to achieve 
nearly Heisenberg-limited metrology even for these sys- 
tems. 

Specifically, we will consider magnetic sensing using 
a single Nitrogen Vacancy (NV) center in diamond 0], 
surrounded by dark spins associated with Nitrogen elec- 
tronic impurities EH • We focus on NV centers since 
their electronic spins (S=l) can be efficiently initialized 
into the Sz = state by optical pumping and measured 
via state selective fluorescence. By applying an exter- 
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nal (static) magnetic field that splits the degeneracy be- 
tween 5*2 = ±1 states and working on resonance with the 
(0 o 1) transition, the NV center can be reduced to an 
effective two-level system Then, the system compris- 
ing one NV center and several N spins is well described 
by Hamiltonian ([T|). 

In Fig. |5]we introduce a control sequence that yields 
an effect equivalent to the circuit in Fig. l.b. The action 
of the pulse sequence can be best understood using the 
well known equivalence between Rarnsey spectroscopy 
and Mach-Zehnder interferometry 18|, ll9[ , where the in- 
terferometer arms describe the central spin state. It is 
sufficient to consider the evolution of each arm sepa- 
rately, replacing Sz by its eigenvalues = {0, 1} and 
describing the evolution in the interaction frame defined 
by the control pulses T?]. Hamiltonian (1) becomes time- 
dependent, with dark spins alternating between /* and 
as shown in Fig. [51 Then, for different halves of the spin 
echo sequence, the coupling Hamiltonian in each arm is 
zero (rris = 0) while for the other halves it has identical 
forms. In the absence of a magnetic field the evolution 
is thus the same along each arm of the interferometer. 
Adding an external field creates a phase shift between 
the two arms. For small field strengths we can then eval- 
uate the phase difference acquired between the two arms 
as a perturbation. For a finite polarization P of the dark 
spins we find 



l + 2P^^sin2W4) 



(4) 



with 9^ = 



dt h{t) - dt b{t) , 6id = f dt b {t). 

Compared to spin-echo based magnetometry the sig- 
nal is increased by the factor in the square bracket while 
the measurement noise is the same, as we still read out 
one spin only. Note that all the dark spins contribute pos- 
itively. For values of the couplings such that |AiT| > tt, 
or strongly coupled environment spins, we obtain a con- 
tribution oc 2nsc (sin^(^i^)) « rise- Each of the n^c 
strongly coupled spins thus gives a contribution of order 
one, irrespective of the sign or exact value of the cou- 
pling. Weakly coupled dark spins (A^r < 1) contribute 
instead with a factor (Xir)'^ /8 and we obtain a total phase 



, where the primed 



1 + |P (n,, -f i (A,t) 

sum is on the weakly coupled spins. In general the sen- 
sitivity enhancement scales as ~ Prise- We thus achieve 
nearly Heisenberg-limited sensing of the external field. 

We next take into account decoherence resulting from 
the interaction with the environment spins as well as de- 
coherence of the dark spins themselves. To evaluate these 
effects, we compare the sensitivity achievable with the 
proposed pulse sequence to that obtained with a spin- 
echo sequence, considering the same system (a sensor 
spin surrounded by the same spin bath) and including 
the effects of decoherence (external perturbations and 
couplings between the dark spins). Once the central spin 




FIG. 2: Method, a) Mach-Zehnder interferometer, showing 
the central spin state and the effective Hamiltonian of the 
dark spins along each arm. b) Pulse sequence for environment 
assisted magnetometry. 



looses phase coherence due to interactions with the bath, 
it is no longer possible to use it for magnetometry. This 
limits the sensing time and consequently the magnetome- 
ter sensitivity. Spin echo (as well as more sophisticated 
decoupling techniques [H, [13, HI, [s^l) can be used to 
prolong the phase coherence of the central spin. Under 
realistic assumptions the coherence time for the pulse se- 
quence presented in this Letter is on the same order of 
the sensor coherence time under spin-echo Tl'. In essence, 
the decoherence rate of an entangled state of the form ([2]) 
is dictated by the internal evolution of the strongly cou- 
pled spins, and the relevant decoherence time is thus the 
time it takes before any of these spins have decohered. 
The same is, however, the case for spin echo sequences: 
if a central spin is strongly coupled to rise dark spins, a 
spin flip of a single dark spin will lead to different evo- 
lutions in the two halves of the spin echo sequence, thus 
decohering the state of the central spin. As the source of 
both decoherence processes is the dipole-dipole interac- 
tion among dark spins, the coherence times of spin-echo 
and our procedure are on the same order and the signal 
amplification obtained with our strategy (Eq. \^ thus 
yields a sensitivity enhancement [33 |. 

To be more quantitative, we analyze the evolution due 
to dipole-dipole couplings in the bath, described by the 
Hamiltonian 



H=\l){l\J2XJl+J2 {sim (5) 



A short time expansion shows that the state fidelity 
at the end of the pulse sequence is given by 1 — 



i<j 



[I'P^) 4(A,-A,)V..., and 1 



128 Z^'^ij ^ 



(7(Af + X]) + 12A,Aj + P2 (A,, + 2Aj) (2A, + X,))+..., for 
the spin-echo and modified pulse sequences, respectively. 
For experimentally relevant values of polarizations, sim- 
ilar reductions in fidelity occur for both pulse sequences 
35[. To go beyond the short time expansions we have 
simulated the signal decay for both spin-echo and the 
proposed pulse sequence (see Fig. [3]). We compared 
these results to the signal decay when no control se- 
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nc=12 




ordered couplings, partial control and decoherence. 
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FIG. 3: Simulations. Normalized signal decay for the pro- 
posed (red, solid lines), spin-echo sequence (blue, dashed 
lines) and no control (black, dotted lines). We assumed per- 
fect delta function pulses. A leading order cluster expansion 
was used [2^. 20 dark spins were randomly placed in a cube 
of side-length \/20 with the sensor spin at the center. We 
set g/iB = 1 [m]'^^'^ [s 
WAHUHA sequences 
echo interval were simulated. Each curve is an average over 
10 Monte Carlo simulations. For simplicity we set P = 0. 



1/2 

' to use dimensionless quantities, 
with He— 8, 12, 25 and 50 cycles per 



quence is applied (for the same environment the decay 
is now described by the dephasing time T2*). We sim- 
ulated a spin bath composed of spin 1/2 paramagnetic 
impurities, undergoing a WAHUHA sequence (which is 
designed to refocus the dipole-dipole coupling of the en- 
vironment spins, but does not cancel out the coupling to 
the external field [B,!!!). 



A different limitation on the sensing time r is set by the 
fact that the orientation of the dark spin will not be static 
as assumed. Dipole-dipole couplings among dark spins 
during each r/4 period of free evolution rotate each spin 
away from the initial direction. This rotation means that 
the spins cease to build up a phase difference between the 
two arms of the interferometer for time scales comparable 
to the correlation time of the dark spin bath t^^. The 
optimum sensing time is thus T ~ min {T|, r^}. Since 
in most systems > T| [17|], the optimum sensing time 



of this pulse sequence is comparable to that of spin-echo 
based magnetometry, thus the sensitivity enhancement is 
roughly the same as the signal strength enhancement. 

In conclusion, we proposed a scheme to enhance pre- 
cision measurement by exploiting the possibility to co- 
herently control the ancillary qubits. In solid state 
implementations we are able to exploit dark spins in 
the bath while preserving roughly the same coherence 
times as in spin-echo based magnetometry. Thus sig- 
nal enhancement leads directly to sensitivity enhance- 
ment. For trapped ion implementations we can use im- 
perfect phase gates and still achieve Heisenberg-limited 
sensitivity. Our method has the potential to be applied 
more generally, using different systems and more sophis- 
ticated pulse sequences [2^. It opens the possibility to 
use a broad class of partially entangled states to achieve 
Heisenberg limited metrology, even in the presence of dis- 
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